Specific heat and thermal conductivity of ferromagnetic magnons in Yttrium Iron 
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The specific heat and thermal conductivity of the insulating ferrimagnet Y3Fe50i2 (Yttrium Iron 
Garnet, YIG) single crystal were measured down to 50 mK. The ferromagnetic magnon specific heat 
Cm shows a characteristic T^'^ dependence down to 0.77 K. Below 0.77 K, a downward deviation 
is observed, which is attributed to the magnetic dipole-dipole interaction with typical magnitude of 
10~* eV. The ferromagnetic magnon thermal conductivity Km does not show the characteristic 
dependence below 0.8 K. To fit the Km data, both magnetic defect scattering effect and dipole-dipole 
interaction are taken into account. These results complete our understanding of the thermodynamic 
and thermal transport properties of the low-lying ferromagnetic magnons. 

PACS numbers: 75.30.Ds, 75.50.Gg 



I. INTRODUCTION 

Recently, the ferrimagntic insulator Y3Fe50i2 (Yt- 
trium Iron Garnet, YIG) draws great attention due to the 
long-range transport ability of spin current.— In these 
experiments, electronic signal can be transferred by spin 
angular momentum of the spin waves in insulating YIG, 
via the spin Hall and inverse spin Hall effects, which pro- 
vides a new method to transfer information by pure spin 
waves. ^'^ In this context, it is important to understand 
the thermodynamic and transport properties of the spin 
waves. 

In spin wave theory for magnets, the quanta of spin 
waves are magnons. There are antiferromagnetic (AFM) 
and ferromagnetic (FM) magnons, which have totally dif- 
ferent dispersion relations, and distinct thermodynamic 
and transport properties. The AFM magnons have lin- 
ear dispersion relation, so both their specific heat and 
boundary-limited thermal conductivity at low tempera- 
ture obey the dependence. For FM magnons, the 
situation is more complex due to the existence of mag- 
netic dipole-dipole interaction (MDDI)<^i^ 

At not very low temperature, the dispersion relation 
E = Dk^ is a good approximation, and the specific 
heat and boundary-limited thermal conductivity of FM 
magnons show the characteristic T^'^ and depen- 
dence, respectively.^"^ The equations are 



15C(5/2)A||5t1-5 
327ri-5i:)i-5 



and 



n^hD 



(1) 



(2) 



where C is the Riemann function, L is the boundary- 
limited mean free pathi^i^Ti^ However, at sub-Kelvin tem- 
perature range, the MDDI 
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with typical order of 10"'* eV, has to be considered for 
FM magnons.*'' It significantly modifies the dispersion 
relation of FM magnons below 1 K, and makes the ap- 
proximate form E — Dk^ no longer valid.** 

The MDDI is long range and anisotropic. It is a basic 
interaction in magnets and critical for many phenomena 
such as the demagnetization factors and the formation 
of domain walls in ferromagnets, the spin ice behavior 
in Ising pyrochlore magnets, and the spin anisotropy in 
ferromagntic filmsJ^"— Theoretically, when the effect of 
MDDI was taken into account for FM magnons, both 
the T*-^ dependence of Cm and dependence of 
changedi^i^ i*"'*^ Yet so far experimental verification of 
this effect on FM magnons is still lacking. 

YIG is an archetypical ferrimagnetic insulator. Its low- 
energy magnetic excitations are FM magnons because at 
low temperature the important spin- wave branch in a fer- 
rite has the same form as in a ferromagnet Ji£i7 ^j^g study 
of this interesting and complex compound started from 
late 1950s, and it has become indispensable for investi- 
gating the properties of magnets since thenJ^ In fact, 
YIG is the first material in which thermal conductiv- 
ity contributed by magnetic excitations was observedJ^ 
At low temperature, its magnon specific heat and ther- 
mal conductivity can even exceed those contributed by 
phononsji^"— Therefore, YIG is an ideal compound to 
study the thermodynamic and transport properties of 
FM magnons. 

Previous lowest temperature for specific heat measure- 
ment on YIG was 1 K i^*'^^ The data of magnon specific 
heat Cm (T) between 1 and 4 K showed the characteristic 
T*-^ dependence The thermal conductivity of YIG 
was roughly measured down to 0.23 The data of 

magnon thermal conductivity Km(T) between 0.23 and 1 
K did not obey the characteristic dependence, and it 
was explained by considering the effect of magnetic defect 
scattering;^i^i2i 

In this paper, we present the specific heat and thermal 
conductivity measurements of YIG single crystal down 
to 50 mK. The magnon specific heat data deviate down- 
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FIG. 1: (Color online) X-ray diffraction pattern for the (332) 
plane of YIG single crystal. Inset: rocking curve of the (664) 

reflection. The two peaks are from the Cu K^i and Ka2 I ■ ' • ' 

radiations, respectively ^ ^ 

-|-1.5 ^1^1.5^ 



ward from the T^-^ dependence below 0.77 K, which is 
attributed to the MDDI effect. Below 0.2 K, the magnon 
thermal conductivity data cannot be fitted by only con- 
sidering the boundary and magnetic defect scatterings, 
and the MDDI effect has to be taken into account. To 
our knowlege, this is the first experimental observation 
of MDDI effect on FM magnon specific heat and ther- 
mal conductivity, giving complete understanding of the 
thermodynamic and thermal transport properties of the 
low-lying FM magnons. 



II. EXPERIMENT 

The YIG single crystal was grown by the optical float- 
ing zone furnace. ^^"^^ The single crystal grew along the 
[332] crystallographic direction, characterized by the X- 
ray diffraction. Ultra-low temperature specific heat mea- 
surement was carried out on a sample with mass 45.55 
mg in a small dilution refrigerator adapted into a Physical 
Property Measurement System (PPMS, Quantum Design 
company). 

Sample 1 (SI) was cut from the single crystal for ther- 
mal conductivity measurements. It is rectangle shaped 
with dimensions 2.08x0.84 mm^ in the plane and 0.63 
mm thick along the [332] direction (the sample growth 
direction). Sample 2 (S2) was obtained by thinning SI 
to 0.23 mm. For both samples, the heat current was 
along the [113] direction. 

Ultra-low temperature thermal conductivity was mea- 
sured in a dilution refrigerator (Oxford Instruments), us- 
ing a standard four-wire steady-state method with two 
Ru02 chip thermometers, calibrated in situ against a 
reference Ru02 thermometer. Four contacts were made 
on the sample surface by silver epoxy. Magnetic fields 
were applied parallel to the heat current. 



FIG. 2: (Color online) Specific heat of YIG single crystal. The 
data between 0.77 and 2.5 K can be fitted by the solid line 
C = 6.7r^-^ + 2.3r^. Below 0.77 K the curve deviates from 
the solid line, which is attributed to the effect of magnetic 
dipole-dipole interaction. Note that the upturn below 0.38 K 
is the Schottky anomaly from nuclear moments. 

III. RESULTS AND DISCUSSION 

The quality of our YIG single crystal was characterized 
by X-ray diffraction (XRD), as shown in Fig. 1. The 
main panel is the XRD pattern of the (332) plane. The 
inset shows the rocking scan curve of the (664) reflection, 
with two peaks from the Cu K^i and Ka_2 radiations, 
respectively. The full width at half maximum (FWHM) 
of the peak from Cu Kai is only 0.07 °, indicating high 
quality of the crystal. 

The specific heat of YIG single crystal below 2.5 K is 
shown in Fig. 2. In the figure, C|T^^^ was poltted as a 
function of T^'^ in order to sperate phonon and magnon 
contributions. Between 0.77 and 2.5 K, the total specific 
heat can be fitted by C = aT^-^^+bT^, in which the first 
term is from the FM magnons and the second term is 
from the phonons. From the fitted coefficient a — 6.7 
/iJ/cm'^ K^-^, we get D — 5.2x10^'^^ J cm^ according 
to Eq. (1). This value is very close to Edmonds and 
Petersen's result D = 5.1x10"^^ J cm^.— 

Below 0.77 K, however, there is an apparent deviation 
from the solid fitting line. As we have described, the 
MDDI will affect the specific heat of FM magnons below 
1 K, by modifying the dispersion relation. The modified 
dispersion relation was proposed as the following forroii 



E{k) = ^{Dk^ - N^huj^){Dk^ - N.huJm + hiOmSm^ Ok) 

(3) 

where 9k is the angle between the magnon wave vector 
and the magnetization direction, N ^ is the z demagne- 
tization factor, and huJm = g^iB^T^M. The dispersion 
relation of Eq. (3) is plotted in Fig. 3. The approxi- 
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FIG. 3: (Color online) The dispersion relation of 
FM magnons at low energy when considering the 
magnetic dipole-dipole interaction effect*^ E{k) = 
\/{Dk^ - N2hbJm)iDk^ - N^hujm + hujm sin^ 6^^). It strongly 
deviates from the E = Dfc^ approximation below the energy 
hum = gfiB'lnM. For YIG, hw^ = 0.32fcfl. 



mate dispersion relation E = Dk^ is only valid for ksT 
> hujm- For YIG, 47rM = 2449 Gs^>^ g = 2,^2 so 
hjj„, = gfiB^T^M = 0.32kB- Therefore, the T^-^ depen- 
dence of Cm should only hold for T > 0.32 K. From 
our experimental data in Fig. 2, Cm shows T^-^ depen- 
dence above 0.77 K. The downward deviation below 0.77 
K should come from the effect of MDDI. Note that the 
upturn below 0.38 K is the Schottky anomaly from nu- 
clear moments. 

Next we discuss the thermal conductivity results. Fig- 
ure 4(a) shows the ultra- low temperature thermal con- 
ductivity of YIG in magnetic fields up to 11 T, plotted 
as k/T vs T. One can see that k is strongly suppressed 
by field, as previously reported »i2ri2S In the insulating 
YIG, the total thermal conductivity can be expressed as 
K — Kpii -|- Km, in which Kp^ and the phonon 

and magnon thermal conductivity, respectively. Since 
Kph is usually not affected by magnetic field, the rapid 
suppression of k with field in Fig. 4(a) should come from 
the reduction of k,„. As we know, the external magnetic 
field H opens a gap A = g/i^TJ in the magnon spectrum. 
When external field is high enough to satisfy gusH ^ 
kbT, there would be no magnon contribution. The sat- 
urated thermal conductivity from iJ = 4 to 11 T shown 
in Fig. 4(a) suggests that there is only phonon contribu- 
tion left below 0.8Kini/>4T, consistent with previous 
reports . ^^1^° 

Therefore the FM magnon thermal conductivity Km in 
zero field can be extracted by subtracting k(4T) from 
k(OT). In Fig. 4(b), Km is plotted as k^/T vs T. For 
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FIG. 4: (Color online) (a) Thermal conductivity of YIG single 
crystal in magnetic fields up to 11 T. In _ff > 4 T, «: tends to 
saturate, (b) Zero field thermal conductivity of FM magnons 
obtained by subtracting k(4T) from «;(0T). The dashed line is 
the fitting of the data between 0.2 and 0.8 K, by considering 
the boundary scattering and magnetic defect scattering. The 
deviation below 0.2 K is attributed to the MDDI effect. The 
solid is the fitting of the data below 0.12 K by including the 
MDDI effect. 



AFM magnons, ballistic boundary-limited Km = oT^ was 
observed below 0.5 K.^ Apparently, for FM magnons in 
Fig. 4(b), there is no such a simple power-law tempera- 
ture dependence. 

Theoretically, the thermal conductivity of magnons is 
calculated by the equation 



fcf 



(A) 



;{VkE)l(fk, (4) 



where i is the mean free path of magnons. At not very 
low temperature, if the approximate dispersion relation 
E = Dk^ of FM magnons is taken and only boundary 
scattering is considered, we get the dependence of Km 
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FIG. 5: (Color online) Thermal conductivity of YIG samples 
SI and S2 in zero field and _ff = 4 T. S2 is obtained by thinning 
SI from 0.63 to 0.23 mm. 



in Eq. (2). Therefore, the anomalous temperature depen- 
dence of Km in Fig. 4(b) must come from the MDDI ef- 
fect or some other scattering mechanism. Previously, ad- 
ditional magnetic defect scattering was considered, and 
the data of Km between 0.23 and 1 K can be well fittedj^S 
In this case, can be expressed as the sum of two terms 



= L- 



(5) 



where L ^ is from boundary scattering and ij^^ is from 
magnetic defect scattering with — ak^. In this way, 
the magnon thermal conductivity is 



^{T) - BT^ 



where a; = E/keT, B 



1 + PT^x^ 



(6) 
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^270^, p - T5^- By using 
this model, our Km data between 0.2 and 0.8 K can also 
be well fitted, with the parameters B = 0.056 mW/K^ 
cm and (3 = 0.15 K~^. Using the obtained value of B, to- 
gether with D = 5.2x10"^^ J cm^, we get the boundary- 
limited mean free path L = 32.7 /im, which we will dis- 
cuss later. 

However, the magnetic defect scattering model can not 
fit the Km data below 0.2 K, as seen in Fig. 4(b). We 
further consider the effect of MDDI. Since the dispersion 
relation Eq. (3) is too complex to do calculation, Orten- 
berger and Sparks chose a simplified dispersion relation^ 



Dk^ + ck 



B- 



(7) 



Substituting this dispersion relation into Eq. (4) results 
in 



Km = BT 



°° x'^csch'^{\x){Tx + c) 
1 + P{Tx -t- c)2 



dx. (8) 



With the B and /3 value obtained above, the Km data can 
be fitted well below 0.12 K, as shown in Fig. 4(b). We 
want to emphasize that we tried to fit the Km data below 
0.8 K with only boundary scattering and MDDI effect, 
without considering the magnetic defect scattering, but it 
did not work. Therefore, the magnetic defect scattering 
mechanism is necessary for explain the low-temperature 
Km data in YIG. Comparing the results of Cm and Km, 
the MDDI starts to affect thermal conductivity at a lower 
temperature than specific heat. The reason may relate to 
the involvement of magnetic defect scattering in thermal 
conductivity. 

Finally, we discuss the boundary-limited mean free 
path L of FM magnons in YIG. From Fig. 4(b), L = 
32.7 /im is obtained by the fitting. This value is one or- 
der smaller than the expected L = 2^J A/n = 727 /im for 
sample SI, with A being the cross section area.— Such a 
phenomenon has been observed previouslyj^^^i^ Fried- 
berg and Harris speculated that it is due to the inner 
boundary of thin layers rich in Fe^'^ inside the samplei^ 
To test this idea, we did the same thermal conductiv- 
ity measurements on sample S2, obtained by thinning SI 
from 0.63 to 0.23 mm. In Fig. 5, the and 4 T data 
of S2 are almost identical to those of SI, indicating the 
same Km in SI and S2. This result shows that Km indeed 
does not change with the sample boundary, therefore the 
actual boundaries are inside the sample, likely the thin 
layers proposed by Friedberg and Harris. These thin 
layers may form during the growing process. 



IV. SUMMARY 

In summary, by extending the measurements of the 
specific heat and thermal conductivity down to 50 mK, 
we investigate the thermodynamic and transport proper- 
ties of the low-lying ferromagnetic magnons in YIG single 
crystal. The deviation of magnon specific heat C„i(T) 
from the characteristic T^-^ dependence below 0.77 K 
is attributed to the effect of magnetic dipole-dipole in- 
teraction. The magnon thermal conductivity Km{T) is 
extracted by subtracting Km(4T) from Km(OT). Below 
0.8 K, Km{T) does not obey the characteristic depen- 
dence due to the magnetic defect scattering. With fur- 
ther decreasing temperature, the magnetic dipole-dipole 
interaction also shows effect on Km{T) below 0.2 K. Our 
work provides a complete understanding of the thermo- 
dynamic and transport properties of the low-lying ferro- 
magnetic magnons. 
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